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Abstract 

Tuberculosis (TB) persists as a public health problem in Iran. Characterization 
of Mycobacterium tuberculosis isolates circulating in this area will contribute to 
understand and control the spread of the strains. The aims of this study were 
to understand the genetic diversity and drug susceptibility of M. tuberculosis 
isolates circulating in Iran and to analyze the relationship between genotype 
and drug resistance. A total of 291 M. tuberculosis isolates collected from TB 
patients were genotyped by spoligotyping. Drug susceptibility testing was per- 
formed using proportion method. Spoligotyping resulted in 75 distinct patterns. 
86.2% of isolates were grouped in 35 clusters while the remaining isolates were 
unique. Ural was found to be the most predominant lineage (34.3%) followed 
by Central Asian strain (CAS) (24%), T (18.2%), Manu2 (7.5%) and Latin 
American-Mediterranean (LAM) (6.1%). The five largest clusters were Ural/ 
Spoligotype International Type (SIT) 127 (15.8%), CAS1/SIT26 (9.2%), Tl/ 
SIT53 (6.1%), T1/SIT284 (5.4%), and CAS1/SIT25 (4.4%). About 5% of isolates 
had multidrug resistance (MDR) and 10% had other resistance. MDR was sig- 
nificantly associated with Beijing strains, but not with Ural family. This study 
highlights dominance of Ural, CAS, and T families in Iran. Biogeographic speci- 
ficity of CAS and T families to border provinces of Iran including Sistan-Balu- 
chestan and Kermanshah, respectively, suggested that this family strains might 
be transmitted fi-om these regions to other provinces of the country. 



introduction 

Tuberculosis (TB) continues to be an important public 
health problem mainly in developing countries. It is esti- 
mated that one-third of the world's population has been 
infected by Mycobacterium tuberculosis. Each year, there 
are more than nine million new TB cases and almost two 
million deaths worldwide despite the fact that it is curable 
with early detection and prompt treatment (WHO 2007). 
Iran has a population of 75 million inhabitants and shares 



geographic borders with three countries where TB is 
endemic, that is Pakistan, Iraq, and Afghanistan. 
Although the reported TB incidence rate for Iran has 
decreased from 36/100,000 in 1990 to 17/100,000 in 2010 
(WHO 2011), TB control remains a priority among pub- 
lic health policies. Molecular typing of M. tuberculosis 
strains is important to identify dominant strains associ- 
ated with outbreak and drug resistance and also to trace 
the transmission chains (Bazira et al. 2011). In general, 
M. tuberculosis isolates with unique fingerprints are most 
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commonly associated with endogenous reactivation 
whereas, clusters of isolates with identical deoxyribonu- 
cleic acid (DNA) patterns have been attributed to recent 
infection (Kaeto-Maeda and Small 2000; Viana-Niero 
et al. 2001). Restriction fragment length polymorphism 
(RFLP) analysis with 7S6110 probe is considered as gold 
standard method for genotyping of M. tuberculosis strains 
(Jiao et al. 2008). However, this method is time-consum- 
ing and labor intensive. It is also unreliable for typing of 
strains with the low copy numbers (fewer than six) of 
7S6110 (Christianson et al. 2010). Spoligotyping is based 
on DNA polymorphisms within the direct repeat (DR) 
locus of M. tuberculosis complex. The DR locus consists 
of a series of directly repeated sequences of 36 bp, inter- 
spersed with 35-41 bp nonrepetitive spacer sequences. 
This method involves a polymerase chain reaction (PGR) 
amplification of the DR locus and hybridization to a 
membrane containing a series of DNA probes represent- 
ing each of the unique spacer sequences in the DR locus. 
It is a simple and rapid method requiring much less DNA 
quantities than 7S6110-RFLP (Kamerbeek et al. 1997). 
Furthermore, recent transfer of spoligotyping technique 
from membrane to a microbead format has even brought 
about more reliable results because of the better detection 
of some spacers and lower error rate in transcription and 
interpretation of the results. Microbead-based spoligotyp- 
ing has enabled analyzing larger number of isolates per 
run which makes it suitable for high-throughput 



genotyping (Cowan et al. 2004; Zhang et al. 2010; Abadia 
et al. 2011). 

Emergence and spread of multidrug resistant (MDR) 
strains of M. tuberculosis represent the greatest threat to 
TB control programs. The design of approaches for the 
management of MDR-TB depends on awareness of devel- 
opment and subsequent spread of drug-resistant isolates 
(Mistry et al. 2002). It is indicated that MDR-TB out- 
breaks have been associated with Beijing and Haarlem 
genotypes (Glynn et al. 2002; Mardassi et al. 2005). 

The specific objectives of this study were to investigate 
population structure of M. tuberculosis isolates circulating 
in five provinces of Iran over a period of 2 years (2010- 
2012) by spoligotyping, and to evaluate the rate of drug 
resistance. In addition, possible association between strain 
type and drug resistance was also studied. 

Material and Methods 

Patient population and bacterial isolates 

This study was conducted at the Department of Microbiol- 
ogy, Tehran University of Medical Sciences, Tehran, Iran. 
A total of 291 M. tuberculosis strains cultured from the 
same number of TB patients were included in this study. 
The isolates were obtained on a survey between November 
2010 to July 2012 from residents of the five different prov- 
inces of the country including Tehran (n = 110), Alborz 
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{n = 14) Sistan-Baluchestan, (n = 89), Hormozgan 
(n = 46), and Kermanshah (n = 32). Location of these 
provinces is shown on a map of Iran (Fig. 1). In total, 270 
isolates were from pulmonary sites (Sputum [n = 165] 
and Bronchoalveolar lavage [n = 105]) while 21 isolates 
were from extrapulmonary sites (Pleural fluid [n = 5], 
wound [n = 8], urine [n = 5], and cerebrospinal fluid 
[n = 3]). Sample processing was performed based on spec- 
imen type. Nonsterile samples were digested and decon- 
taminated with N-acetyl-L-cycteine sodium hydroxide 
(NALC-NaOH) method (Kubica et al. 1963; Kent and 
Kuniba 1985). In brief, 10-15 mL of specimen was mixed 
with equal volume of NALC-NaOH solution and homoge- 
nized with a vortex mixer for 15-20 sec. After incubation 
for 15 min, the mixture was neutralized with phosphate 
buffer, pH 6.8 and concentrated by centrifugation at 3000^ 
for 15 min. The sediment was resuspended in 0.5 mL of 
sterile phosphate buffer, pH 6.8 and cultured on Lowen- 
stein-Jensen (LJ) slant. Sterile specimens were processed 
and cultured without prior decontamination. LJ slants 
were incubated for up to 8 weeks. M. tuberculosis isolates 
were identified using standard biochemical tests, including 
catalase activity, production of niacin, nitrate reduction, 
pigment production, and growth rate (Kubica 1973). 

Drug susceptibility testing 

Drug susceptibility testing against first-line drugs was per- 
formed according to proportion method on LJ media with 
the following drug concentrations: isoniazid (INH), 
0.2 /ig/mL; rifampicin (RF), 40 /.(g/mL; streptomycin (SM) 
4 /(g/mL, and ethambutol (ETB) 2.0 /xg/mL (Canetti et al. 
1963, 1969). Briefly, M. tuberculosis colonies were placed in 
thick-walled screw-capped glass tube containing 10 sterile 
glass beads (3 mm in diameter). Tubes were shaken for 
20-30 sec. Distilled water was added slowly untU the opac- 
ity of the bacterial suspension matched with that of 1 mg/ 
mL of Mycobacterium bovis Bacillus Calmette-Guerin 
(BCG). Ten-fold serial dilutions were prepared down to a 
dilution of 10"^ and 0.1 mL of dilutions 10"^ and 10"^ 
was inoculated on to the control (drug-free) and drug-con- 
taining media. The extent of growth in the absence or pres- 
ence of drug was compared and expressed as a percentage. 
Strains were considered as sensitive when the growth rate 
was <1% in comparison to control. MDR was defined as 
resistance to at least INH and RF (WHO 2010). This pro- 
ject has been approved by Tehran University of Medical 
Sciences Ethical Committee (Project No. 15356). 

DNA extraction and spoligotyping 

Deoxyribonucleic acid was extracted from mycobacterial 
colonies growing on LJ medium using the method of 



Wilson (1987). Spoligotyping was performed as previously 
described by Kamerbeek et al. (1997) using membrane 
provided from National Tuberculosis Reference Labora- 
tory, Bilthoven, the Netherlands. In brief, DR region of 
mycobacterial genome was amplified using PCR with 
primers DRa 5' - GGT TTT GGG TCT GAC GAC -3' 
(biotinylated at 5' end) and DRb 5'-CCG AGA GGG 
GAC GGA AAC-3' (Metabion, Martinsried, Germany). 
The PCR products were subsequently hybridized to a set 
of 43 oligonucleotides covalently bound to a membrane. 
The hybridized PCR products were incubated with 
streptavidin-peroxidase conjugate (Roche, Mannheim, 
Germany) and then detected by chemiluminescence sys- 
tem (Amersham ECL detection kit; GE Healthcare Lim- 
ited, Little Chalfont, Buckinghamshire, U.K). The X-ray 
film was developed using standard photochemical proce- 
dures after 1-30 min exposure. Two membranes were 
purchased for genotyping of 291 samples. AH samples 
were genotyped in duplicate by using two different mem- 
branes. For repeated use of membranes each membrane 
was regenerated eight to nine times by two times incuba- 
tion in 1% sodiumdodecyl sulphate (SDS) at 80°C for 
30 min, and subsequent washing in 20 mmol/L ethylene- 
diaminetetraacetic acid for 15 min at room temperature 
(Kamerbeek et al. 1997). DNA extracts of M. tuberculosis 
H37Rv and M. hovis BCG were used as positive controls. 

Spoligotyping data analysis 

We compared our spoligotyping results in binary format 
with the SITVIT2 database (Pasteur Institute of Guade- 
loupe) an updated version of the international spoligotyp- 
ing database, SpolDB4.0 (Brudey et al. 2006) (available 
online at http://www.pasteur-guadeloupe.fr:8081/SITVIT- 
Demo). A spoligo type-based dendrogram was generated 
by the unweighted pair group method with arithmetic 
averages (UPGMA) calculation. A cluster was defined as 
two or more isolates from different patients with identical 
spoligotype patterns. We defined unique (nonclustered) 
spoligotypes as those which did not cluster with any other 
sample in our study. Orphan was defined as a unique 
spoligotype pattern not described in the SITVIT2 data- 
base. Statistical analysis was performed using chi-square 
test to determine statistical associations between strain 
types and specific parameters and values of P < 0.05 were 
considered statistically significant. 

Results 

Spoligotyping analysis 

Spoligotyping of 291 isolates produced 75 distinct 
spoligotype patterns. Overall, 251 isolates (86.2%) were 
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Table 1. Spoligoptyes distribution 
groups. 


in clustered 


and nonclustered 


Spoligotypes 


Patterns 


No. of isolates (%) 


Previously identified spoligotypes 


54 


252 (86.5) 


Clusters 


31 


229 (78.6) 


Unique 


23 


23 (7.9) 


Orphan spoligotypes 


21 


39 (13.4) 


Clusters 


4 


22 (7.5) 


Unique 


17 


17 (5.8) 


Total clusters 


35 


251 (86.2) 


Total unique 


40 


40 (13.7) 


Total 


75 


291 



classified into 35 clusters containing 2-46 isolates while 
the remaining 40 isolates (13.7%) did not form clusters 
(unique patterns) (Table 1). Spoligotyping grouped 278 
(95.5%) isolates in 10 families including Ural (former H4 
sub-lineage, n = 100, 34.3%), Central Asian strain (CAS) 
(CASl and CAS2, n = 70, 24%), T (Tl, T2, T3 and T4, 
n = 53, 18.2%), Manu2 (n = 22, 7.5%), Latin American- 
Mediterranean (LAM) (LAMl, LAM2, LAM6, LAM9, 
n = 18, 6.1%), East African Indian (EAI) (n = 4, 1.3%), 
Beijing (n = 3, 1%), U (n = 3, 1%), Haarlem (H3, n = 3, 
1%), and H37Rv (n = 2, 0.6%). Twenty-one spoligotype 
patterns (n = 39 isolates) were defined as orphan, among 
them orphan 1, 3, 5, 6, 7, 8, 14, 19, and 21 (n = 25 iso- 
lates) showed high similarity to Ural genotype and were 
classified with Ural linage. Similarly, orphan spoligotype 
10 (m = 1 isolate) was classified as EAI family based on 
high similarity to this family. The remaining 11 orphan 
spoligotypes (including orphan 2, 4, 9-13, 15-18 and 20, 
M = 13 isolate) were not found in the database SITVIT2 
and did not show similarity to any other known lineages 
and were considered as possibly novel genotypes. Table 2 
depicts the spoligotypes of 10 most prominent 



spoligotypes with their SIT number in SITVIT2 database. 
Table 3 shows the spoligotype description of 39 orphan 
isolates obtained in this study. The most frequent spoligo- 
type in our population was Ural/SIT127 (n = 46, 15.8%) 
followed by CAS1/SIT26 (n = 27, 9.2%), T1/SIT53 
in = 18, 6.1%), T1/SIT284 (n = 16, 5.4%), CAS1/SIT25 
in = 13, 4.4%), and Manu2/SIT54 (n = 12, 4.1%). An 
autoradiograph of spoligotyping representing hybridiza- 
tion pattern of amplified DR locus from M. tuberculosis 
isolates is depicted in Figure 2. The UPGMA tree of 
M. tuberculosis genotypes is illustrated in Figure 3. A 
comparison of genotype distribution among five prov- 
inces indicated that prevalence of Ural family (n = 58/ 
124, 49.7%) was higher in Tehran and Alborze, while 
CAS superfamUy was more frequent (n = 40/89, 45%) in 
Sistan-Baluchestan. T (n = 16/32, 50%) and LAM (12/46, 
26%) clades were prominent in Kermanshah and Hor- 
mozgan, respectively. Overall, three Beijing strains (SITl) 
were identified consisting of about 1% of the isolates 
studied. Two of these isolates were from Tehran province 
and one was from Hormozgan province. 

Patient characteristics 

A total of 291 M. tuberculosis strains were isolated from 
135 women aged between 9 and 88 and 156 men aged 
between 2 and 89. The average age of studied group was 
56 years. Patients in the group aged >50 years had the 
highest number (52%, n = 152) of isolates, followed by 
those aged between 30 and 50 with 30.5% (« = 89) of the 
isolates. 



Drug resistance patterns 

Of the 291 strains in this study 84.8% (n = 247) were 
sensitive to all tested agents. Among the 44 drug-resistant 



Table 2. Spoligotypes of 10 most prevalent genotypes of Iran. 



SIT' 



Family 



No. (%f 



Spoligotype pattern 



127 

26 

53 

284 

25 

54 

New 

754 

64 

20 



Ural 
CASl 
Tl 
Tl 

CASl 
Manu2 

Ural (Orphan7) 
CASl 
LAM6 
LAMl 



46 (15.8) 
27 (9.2) 
18 (6.1) 
16 (5.4) 
13 (4.4) 
12 (4.1) 
12 (4.1) 
11 (3.7) 
8 (2.7) 
5(1.7) 



■■■■■■■□□□■□□□□■■I 
■□□□□□□□□□□□□■■■■I 
■■■■■■■■■■■□□□□■■I 
■■■■■■■■■■■□□□□■■I 

■□□□□□□□□□□□□■■□□I 



■■■□□□□□□□□□□□□■■I 



^SIT from SITVIT2. 

^Representing spoligotype families annotated in SITVIT2. 
^Number of strains. 

^(■) Presence of spacer; (□) absence of spacer. 
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Table 3. Spoligotype patterns of orphan isolates of Mycobacterium tuberculosis in Iran. 



Name of spoligopattern 


No. of isolates (%) 


Spoligotype pattern 




Orphan seril ^ 








Orphani 


1 (0.3) 






Orphans 


2 (0.6) 






Orphans 


5 (1.?) 






Orphans 


1 (0.3) 






Orphan? 


12 (4.1) 






Orphans 


1 (0.3) 






OrphanU 


1 (0.3) 




Orphani 9 


1 (0.3) 






Orphan21 


1 (0.3) 




Orphan seri2^ 








OrphanIO 


1 (0.3) 






Orphan seri3^ 








Orphan2 


1 (0.3) 






Orphan4 


1 (0.3) 






Orphan9 


3 (1) 






rirnhsn 1 1 

Wl Lt[ Idl 1 1 1 


1 {0 3) 




Orphani 2 


1 (0.3) 




Orphani 3 


1 (0.3) 






Orphani 5 


1 (0.3) 






Orphani 6 


1 (0.3) 






Orphani? 


1 (0.3) 






OrphanIS 


1 (0.3) 




Orphan20 


1 (0.3) 







'Orphan spoligotypes related to Ural lineage. 
^Orphan spoligotype related to EAI lineage. 

^Orphan spoligotypes not identified as genotype of any phylogenetic lineage. 
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Figure 2. Spoligotyping autoradiograph 
representing hybridization patterns of amplified 
DNA samples from 39 Mycobacterium 
tuberculosis clinical isolates, 2 M. tuberculosis 
H37Rv and 1 Mycobacterium bovis BCG. 
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UPGMA-Tree, Spoligo; Categorical 




Figure 3. Spoligotype-based dendrogram of Iranian isolates (n = 291) 
of isolates enclosed in parentheses. 



Boxes with yellow background show spoligotypes SIT number and number 
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isolates, 15 (34%) were MDR and remaining 29 isolates 
showed mono drug resistance for INH (n = 4), RF 
in = 2), SM in = 21), and ETB (n = 2). 7.5% (« = 3/40) 
of nonclustered isolates were drug-resistant compared 
with 16.3% (n = 41/251) of clustered isolates. This differ- 
ence was not statistically significant (P-value = 0.061 
[Pearson's chi-square test]). 

Association between genotype and drug 
resistance phenotype 

In an investigation between M. tuberculosis genotypes and 
drug resistance we found that Beijing genotype was highly 
associated with MDR. Actually, all three Beijing strains 
showed MDR phenotype. Only 5% of Ural family isolates 
were MDR in this study. This figure is not significantly 
different from that of other circulating genotypes in our 
setting (P = 0.99, chi-square test) (Table 4). 

Discussion 

In this study, we provided information about the geno- 
typic biodiversity and drug resistance of M. tuberculosis 
strains collected from five provinces of Iran. About 86.2% 
of our isolates were clustered by spoligotyping. Previous 
spoligotyping studies performed in Iran (Farnia et al. 
2006; Ramazanzadeh et al. 2009) showed that Haarlem 
was the most frequent lineage. East African Indian was 
the most predominant family in another study reported 
from Iran (Velayati et al. 2006). However, the majority of 
M. tuberculosis strains presented in this study was charac- 
terized with Ural family. Ural (former H4 sub-lineage) 
lineage has been recently excluded from the Haarlem fam- 
ily based on not harboring mgtC545iCGC -> CAC) 
mutation (Abadia et al. 2010). The CAS family was 
ranked as the second most prevalent spoligotype. It is 
overrepresented by prototype SIT26 and was causative of 
TB in 24% of the cases. Principally, 57% (« = 40/70) of 
this family were from Sistan-Baluchestan which is the 
largest province of Iran with the highest rate of TB in 



the country. This province is located in the southeast of 
the country and shares geographic borders with Afghani- 
stan and Pakistan where the CAS family is dominant. In 
its recent history Iran has been ranked as world's largest 
refugee haven, mainly for Afghans and Iraqis. Iran also 
hosts some refugees of various nationalities, the majority 
of these being Pakistanis, Afghans, Tajiks, Bosnians, and 
Azeris. Also, at the end of 2005, the UN Refugee Agency 
estimated that Iran was host to the third largest refugee 
population in the world, with a total of 716,000 refugees 
(http://www.migrationinformation.org). Because of the 
high prevalence of CAS family (especially SIT26) in Paki- 
stan (Hasan et al. 2006; Tanveer et al. 2008) it is possible 
that this family strains were transported from Pakistan to 
Sistan-Baluchestan. However, the proof of this hypothesis 
is difficult as more population-based epidemiologic stud- 
ies are required. The prevalence of Beijing strains shows 
geographic variation. High rates of Beijing genotype have 
been reported from Asian countries (China, Hong Kong, 
and Thailand >50%) (Glynn et al. 2002), while in Wes- 
tern European and Czech Republic (<6%) and Latin 
America (<1%) its prevalence is low (European Concerted 
Action on New Generation Genetic Markers and Tech- 
niques for the Epidemiology and Control of Tuberculosis 
2006). Previous research from Iran showed that the per- 
centage of Beijing genotype was 20.5% (Farnia et al. 
2006). However, this genotype comprised only 1% of our 
studied strains. The explanation for this difference may 
be that the strains used in the previous study were 
obtained from patients with MDR-TB which is highly 
associated with Beijing genotype. Most M. tuberculosis 
isolates recovered in Iran were susceptible to all tested an- 
tituberculosis drugs. In this study, about 10% of the 
isolates showed mono drug resistance and 5% of isolates 
were MDR. Association between Beijing genotype and 
MDR-TB varies worldwide. While researches from Russia, 
Cuba (Glynn et al. 2002), and Pakistan (Tanveer et al. 
2008) reported such an association, it has not been noted 
in reports from Indonesia, Western Sweden (Brudey et al. 
2004), and Hong Kong (Glynn et al. 2002). Analysis of 



Table 4. Association between bacterial genotype and drug resistance. 



No. of strains (%) 



Characteristic 


Ural 


CAS 


Beijing 


Orphan 


Others' 


All spoligotypes 


INH resistance 


3(3) 


0 


0 


0 


1 (0.9) 


4(1.3) 


RF resistance 


1 (1) 


0 


0 


0 


1 (0.9) 


2 (0.6) 


SM/ETB resistance 


12 (12) 


4 (5.7) 


0 


0 


7 (6.6) 


23 (8) 


MDR 


5(5) 


1 (1.4) 


3 (100) 


0 


6 (5.7) 


15(5) 


Sensitivity to all drugs 


79 (79) 


65 (92.8) 


0 


13 (100) 


90 (85.7) 


247 (84.8) 


Total 


100 


70 


3 


13 


105 


291 



'includes T, Manu2, tAM, EAI, U, Haarlem, and H37Rv families. 
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association between bacterial genotype and drug resis- 
tance showed relationship between Beijing isolates and 
MDR-TB. However, the proportion of this genotype in 
this study population is too small. Therefore, a greater 
study is required to confirm this finding. In spite of the 
high prevalence, only five out of the 100 (5%) Ural type 
isolates were MDR in this study. This prevalence was not 
significantly different from that of the other circulating 
genotypes in our population. The third most predomi- 
nant family in this study was ill-defined T family (in 
which spacers 33-36 are absent), which is also reported 
from western Sweden (Brudey et al. 2004), Brazil (Men- 
des et al. 2011) and Pakistan (Tanveer et al. 2008). The 
strains of this family comprised about 18% of all isolates 
and 50% of isolates from Kermanshah Province. This 
province is a part of region called Iranian Kurdistan, 
which is an unofficial name for the parts of Iran inhab- 
ited by Kurdish people and has borders with Iraq and 
Turkey. In several spoligotyping studies carried out in Tur- 
key (Gencer and Shinnick 2005; Kisa et al. 2007; Aktas 
et al. 2008) ill-defined T family (especially SIT53) was 
reported to be the most frequent genotype. This family was 
demonstrated as the predominant family in a study per- 
formed in Iraq (Merza et al. 2011). This correlation sug- 
gests a possible role for migration, trade and tourism in 
distribution of T family among these three countries ("Tur- 
key, Iran ready to bolster tourism" Turkish daily news. 19 
June 2006; http://www.thenational.ae/business/travel-tour- 
ism/iraqi-citys-religious-tourism-set-to-grow). 

In conclusion, this study gives preliminary information 
about major M. tuberculosis strains circulating in Iran and 
their drug resistance pattern. The main Iranian families 
found in this study were Ural, CAS, and T. All Beijing 
strains were MDR, whereas the most frequent lineage, 
Ural did not show significant association with antituber- 
culosis drug resistance. We suggested that CAS and T 
families may be transmitted from bordering provinces of 
the Iran to other regions of the country. Genotyping with 
other powerful methods and studying other provinces in 
Iran wiU be the next step of this ongoing work to provide 
better understanding of the epidemiology of TB within 
Iran. 
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